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We consider Raman scattering of coupled two leg spin ladders in the Fleury-Loudon regime. We 
derive the dependence of the intensity with polarization of the incoming light and temperature and 
discuss the effect of interaction of the elementary excitations on the shape of the spectrum. We 
show that Raman scattering spectrum is sensitive to the effective dimensionality of the coupled 
ladder system, making Raman scattering useful to discriminate true ladder systems from quasi two 
dimensional spin gap systems. 
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The study of spin ladder compounds! has become a 
subject of intense theoretical and experimental activitycl 
in the recent years. Insulating two leg spin-ladder com- 
pounds present a spin gap state that has been probed ex- 
perimentally by thermodynamic (specific heat and mag- 
netic susceptibility) measurements as well as dynami- 
cal measurements such as neutron scattering or NMR. 
Another technique to study the dynamics of magnetic 
excitations in low dimensional antiferromagnets is Ra- 
man scattering. This technique has been used ia par- 
ticular to probe spin- 172 chainsa, spin 1 chainsu and 
spin Peierls systemsETu. The extended interest to- 
wards this experimental technique is due to its sensi- 
tiveness to singlet excitations, that render it complemen- 
tary to neutron diffraction experiments. There exists 
at present a certain amount of literature on the theory 
of Ramaix. scattering from dipierized spin chains, both 
analyticaoEj and numericaJtiLin relation to experiments 
on Spin-Peierls compoundsETEl. In the case of spin lad- 
der materials, although there are quite a few experimen- 
tal investigational3il3 by Raman scattering, there is a 
scarcity of theoretical results on this subject. While some 
numerical calculations are available to describe the exper- 
imental dataHl, no analytical expressions of the Raman 
intensity for 2n-legged spin-1 /2 ladders existed. In all the 
theoretical work on spin ladders, the analysis of mag- 
netic Raman— scattering is based on the Fleury-Loudon 
HamiltoniarJiia Eil, that describes the interaction of pho- 
tons with magnetic excitations. This limits the investi- 
gations to the non-resonant regimeEH. In a previous re- 
porttSI, we have derived the relevant expressions for two 
leg ladders with weak and strong interladder coupling. 
In our reportta, we found that for weak coupling case, 
the Majorana fermions description of the spin ladderEj 
lead to a cusp in the Raman intensity at twice the gap. 
This was in disagreement with experiments that show 
a peak rather than a cusp in the Raman intensity for 
a frequency twice the gap. On the other hand, ha—the. 
strong coupling case, the Bond Operator TcchniqueEl~c3 
(BOT) predicted correctly the presence of the peak at 



twice the gap. This suggests that the strong coupling de- 
scription is more adequate to describe Raman scattering 
experiments in real spin ladder systems. In our report, 
some important aspects of Raman scattering in spin lad- 
der were not considered. First, we restricted ourselves 
to zero temperature. However, the spin gap is tempera- 
ture dependent and thus temperature should modify the 
Raman intensity. A second aspect that was not consid- 
ered is the effect of interactions between the magnetic 
excitations. Such interactions break the coherence that 
gives rise to the peaks and it is important to understand 
how they modify the ideal spectrum of Ref. A last 
aspect is the influence of iaterladder excitations. There 
are spin ladder materialsEStZl in which the interladder 
coupling is by no means negligible such as KCUCI3 or 
VO2P2O7. This affects strongly the neutron scattering 
intensity compared to the one of a single ladder. This 
two-dimensionality should also affect the Raman spec- 
trum. Understanding how the Raman intensity is af- 
fected by dimensionality is important to discriminate a 
one-dimensional from a two-dimensional compound. A 
possible application could be to the CuHpCl compound 
for which— there exists some evidence of two dimensional 
behavioral (See however Refs. |f-fl|). 

In the present paper we deal with the calculation of the 
Raman scattering cross section of coupled spin ladders. 
After recalling some basic results on the Fleury-Loudon 
theory of magnetic Eaman scattering, and the Bond Op- 
erator TechniqueEjta (BOT) we give the general deriva- 
tion of the Raman operator for coupled ladders in Sec. 
I In Sec. H we give .an expanded discussion of the two 
leg Heisenberg ladderliS. We discuss the temperature de- 
pendence of the Raman intensity of the ladder. We also 
consider the effect of interactions on the peaks in the Ra- 
man spectrum and derive a scaling form of the intensity 
valid in the vicinity of the edges of the Raman spectrum. 
Having illustrated the technique on and the p hysics on 
the simple isolated ladder case, we turn in Sec. [II B 2 to 
the case of coupled ladders. First, we consider a system 
of two coupled ladders, for which the Raman spectrum 
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presents four peaks. The influence of temperature and 
interactions on the spectrum is the same as in the case of 
the isolated ladder. Then, we consider the case of a two 
dimensional array of ladders. We show that the spec- 
trum of this system is markedly different from that of 
the isolated ladder, and discuss the effect of the interac- 
tions of magnetic excitations on the singularities of the 
two-dimensional spectrum. 



I. THE LOUDON-FLEURY THEORY FOR THE 
PERIODIC ARRAY OF LADDERS 

We consider a periodic array of single-rung two-legged 
Heisenberg ladders (shown in Fig.|l|) whose Hamiltonian 
is 

H = (J\\Si j ■ Sj +lj - + J±S i 2 j ■ S i 2 j+i 



(1) 



where J11 = AJj_ is the interaction along the legs of the 
ladder, J± is the rung interaction and J' ± = \'J± is the 
interladder coupling. The index i and j run along the leg 
and rung direction respectively. In the limit A' = the 
Hamiltonian reduces to a sum of independent single-rung 
ladders, while in the limit A = A' = (strong-coupling) 
the Hamiltonian reduces to a sum over independent two 
spin rungs. In this paper we are interested in the regime 
A'<A<1. 

The interaction of light with the antifeHjemagnetic 
fluctuations is described by Loudon-Fleury'al3lH3 photon- 
induced super-exchange operator 



H R = 5^7ii'(Ei • M(E B • 5ii')Si • a v 



(2) 



1,1' 



where Ei(Es) are the incident (scattered) electric field 
vectors of photon, and S\\' is a unit vector connecting the 
lattice sites 1 and 1', at which the spins Si and Si' are 
located. A derivation of (||) starting from the Hubbard 
Hamiltonian can be found inJleL pl[ 

The Raman cross sectiorHEj can be expressed in 
terms of of the retarded Raman response function as: 



LO1LO2 TI2 



2ttcW m 1 



(3) 



uj\ and LU2 are the frequencies of the incoming and scat- 
tered radiation, respectively, lu — lo 2 — u>i, n\ and ni are 
the refractive indices of the material at frequency respec- 
tively to\ and u>2 ■ V is the volume of the crystal and c the 
velocity of light. The retarded linear response function 
X.r(w) is defined as: 



ret 

Xr 



1 f°° 

(w) = ~ J e^+^'Tr {Z- X e- &B [H R (t), H R (0)}} , 

(4) 



where Z = Tre^^and Hr is the Loudon-Fleury Hamil- 
tonian (|J). 

By inserting the resolution of identity in (^) , the Ra- 
man intensity can be written as 



'-a 11^ 
oc > ( 



-0E n 
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|<*„ \H R \ ^>m)\ 2 5{u - {E n - E m )/h), 



(5) 



where | & n (m)) are eigenstates with energies E n ^ m y Such 
formula can be easily interpreted as a Fermi golden rule 
averaged over the Boltzmann weight. To obtain infor- 
mations on two-magnons scattering processes we should 
perform a symmetry analysis of the matrix elements ap- 
pearing in rtq), and discuss selection rules. Since the spin 
ladder Hamiltonian is invariant under translation along 
the legs, SU(2) rotation, and mirror along the the leg di- 
rection, an eigenstate should be characterized by a (lat- 
tice) momentum defined modulo 2ir/a (where a is the lat- 
tice spacing), a spin and its parity under leg exchange. 
The Raman operator defined in (|J) is rotationally and 
translationally invariant, and still invariant under leg ex- 
change. As a result, the selection rules impose that the 
states I $„) and | ^> m ) have the same spin, momentum 
and parity under leg exchange. This implies in particular 
that at T = 0, transitions will only take place to states of 
total momentum zero, spin zero and same parity as the 
ground state. 

To obtain an expression of the Raman Hamiltonian 
more adapted to perform the calculations, we introduce 
a coordinate system [x, y) where x is parallel to the chain 
direction and y is parallel to the rung direction. In this 
coordinate system, we have 1 = ix + jy = Ej = 

E '1 (cos 61$ + sin 9 iy) and Es = Es(cos9sx + sin 9sy). 
We define 7|| = 7(i,j),(i+i,i)> 7-L = 7(i,2j),(i,2j+i) and 7^ = 
7(i,2j'+i),(i,2j+2)' The expression of the Raman operator 
becomes: 

Hr = EiE s cos 9 S cos 6ij\\ Sij ■ S,+ij + EiE s sin 9j sin 9s 



X I 7-L Si f 2j • Si i2 j+1 + 7j. Sj^j'+l • &i,2j 



+ 2 



(6) 



This expression can be simplified using the form of the 
full Hamiltonian. One has: 



7i 



H R = -±E I E s sm9 I sm9 s H 

'J I 



( 7i J ll A 
-EiE s I 7|| cos 9 S cos 9i -p— sin 9i sin 9 S J &ij ' s »+i 



-E]Es sin 9 1 sin 9s 7_l 



y . ■ Si,2j+i (7) 



1,3 



The part of Hr proportional to H gives no contribution 
to the response functions. 
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In the case of a single ladder, J' ± = and the simplified 
form of Hji is: 



, / 7-L^|| 

H R — I711 cosOi cosOs — — sin 81 sin 8 s 

x EjEs 2_j 53 ' 



(8) 



Some remarks on the coefficients j a , (a = ||,_L) ace in 
order here. According to perturbative calculations^ for 
a half filled Hubbard ladder, one should have 



J±/J± =l'±/J± =1\\/ J \\ 



(9) 



With this relation, (g) is recovered when taking the limit 
J' ± — > in (0). In such case, the simplified form of the 
Raman Hamiltonian is: 



H R = 7||-Ej-Es cos(0/ + S ) 2^ Si j • Si+ij 



(10) 



This form has been previously obtained in Ref. |34|. It 
leads to the same Raman intensity for 81 = 9s = a«d 
81 = 83 = 7r/2 in disagreement with experimental dataEa. 
It was shown that treating the 7's as phenomenological 
coefficients and assuming j a oc \/Ja lead to better agree- 
ment with experiments^. In the following, we will treat 
the 7's as phenomenological coefficients and then spe- 
cialize to the simplified expression that result when (||) 
is assumed. 

To analyze the Raman susceptibility, we use the Bond 
Operator Representation (BOT)cJ'E£l of .quantum S=l/2 
spins used by Gopalan, Rice and SigristEa in their mean 
field approach to spin ladders. In this representation, 
one starts from weakly coupled rungs, i.e. A = 0, and 
introduces on each rung a singlet and three triplets 
t' a (a = x,y,z) boson creation operators, that span the 
Hilbert space of a single rung when acting on a vacuum 
state. Since the rung can be in either the singlet or one 
of the triplet states, the condition : 



» + s + J^tttc, = 1 



(11) 



has to be satisfied by the physical states. The representa- 
tion of the spins Si^j and S^j+i in terms.o£ these singlet 
and triplet operators, is derived in Rcf.ESEj. It reads: 



,2j 



(12) 
(13) 



The Hamiltonian Eq. (|l|) can be rewritten in the BOT 
as: 



XJ ± s 2 



^'J-L-T 5Z + 4,3,o) (*».J+l.o + 4,3+l,a) 



i,3=l,2 



1,3 



53 ( — T^i s2 — M^ 2 + M ) + quartic terms. (14) 



The quartic terms in t, v in Eq. (14) describe interactions 
between triplet excitations. Their explicit form will not 
be needed here. 

Substituting the operator representation of spins 
Eq. rtl3) into the Raman operator for coupled ladders 
Eq. (Mj, one ends up with the following expression: 



H R = EiEs I 7|| cos t>s cos t>i =-j-^ sin f 7 sin ( 



i+l,2j,t» + U+l,2j,a ){t 



EiEssm8ism8s (-fx - ^) £ J'Uate ( 15 ) 



».3.a 



In Fourier space, this expression becomes: 
EiEs 



H R = 



k , a 



where: 

^k-(7±- i y i )sin0/sin^ 

■'i. 

7' J\\ 

+ 2s 2 cos(fc x a)(7|| cos 9i cos 9s -, — sin#/ sin6>s) 

•*± 

rJ J,, 

•Bk = s 2 7|| cos 0/ cos 9 S - s 2 sin 81 sin 6> s (17) 

In particular, if 7||/</|| = 7_i_ /■/_!_ = l'±/J±, wc get the 
following relation between the coefficients 

Ak = 2£?k cos(fc x a) = 2s 2 7|| cos k x a cos(#/ + 9s). (18) 

The above expressions will be used in the following of 
the paper to predict the Raman intensity. 



II. THE SINGLE RUNG LADDER 

A. Mean field theory 

We start by considering a single ladder of two strongly 
coupled antiferromagnetic 5 = 1/2 Hcisenberg chains, 
whose Hamiltonian is 

H = J|| 53 (Sl.tSl^+l + S2,iS2,i+l) + J_L 53 Sl,-i ' ^2,i 

i i 

(19) 
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where Jii = \J± > and J± > denotes the intra- and 
inter- chain antiferromagnetic interactions, respectively. 
Substituting the BOT operator representation of spins 
( |l3| ) into the original Hamiltonian, one ends up with an 
Hamiltonian quartic in boson fields. Treating the singlet 
operator in a mean field approximation and neglecting 
interactions among the triplets, one obtains_the following 
Hamiltonian quadratic in triplet operatorso: 



H 



MF 



Js 2 



y^fti.q + **,a)(*i+l,e< + k+l,a)- 



(20) 



The chemical potential term [i guarantees that the condi- 
tion ( p"l| ) is satisfied on average. The parameters /z and s 
are determined in a self-consistent way by the minimiza- 
tion of the free-energy. The self consistent equations to 
be solved are: 



where 



77 dk 
Am 



coth 



2 



V 1 + d cos k 



Vi- 


- dcosk 


4A , 


r dk 




— coth 


~d J 


47T 







yl + d cos k = 0. 



d = 



(k- JL)' 



V f + d cos fc 

(21) 



(22) 



After some simple manipulations, from equations (|2 
we obtain the following equation for d 



d = A 



3- 



coth 



7T Vl + dcosk 



(23) 



When this equation is combined with the first equa- 
tion in (pl|), we can solve them numerically to obtain s 2 
and d. Only in the zero temperature limit these equa- 
tions decouple completely and reduce to the equations 
(2.19) of Ref. 25. For the purpose of this paper, i.e. the 
analysis of the Raman intensity, we discuss the tempera- 
ture dependence of the singlet order parameter s (T) and 
d(T). The numerical analysis shows (see Fig.(^)) that at 
high enough temperature s is very small and increases as 
temperature is lowered. Such behavior is expected since 
at high enough temperature bosons cannot condense. In 
this high temperature regime where s(T) = the mean- 
field Hamiltonian ( p0| ) commutes with the Raman op- 
erator, leading to zero Raman intensity from magnetic 
scattering. For low enough temperature, s(T) ^ and 
a Raman signal appears. As temperature is lowered, the 



bandwidth of triplet excitations increases and the width 
of the magnetic Raman scattering increases. Let us point 
out that the mean field calculation of s(T) is valid only for 
a low density of triplets i. e. for low enough temperature. 
At high temperature, the triplet-triplet interactions can- 
not be neglected. However, the mean-field theory gives 
nevertheless a qualitatively correct picture of the devel- 
opment of Raman intensity. The Raman intensity (see 
Appendix) is given by: 



Iiiixk(w) = C 2 (9 I ,9 S ) / dfccoth 
- 8{lo + 2uj k )) 

where 



2k B T 



C(& I ,9 S ) = - 



7j_ sin t>i sin its 



1\\J± 



COS Vl COS 6*5 



(25) 



It measures, up to a matrix element, the density of 
states of the triplet excitations. The Raman scattering 
spectra thus displays two peaks at its edges, the first one 
at energy u) = 2lo^ = 2A S corresponding to the bottom 
of the triplet band, and the second one at u> = 2ll>q, corre- 
sponding to the top of the triplet-baud . Such peaks have 
been observed experimentallyEjO'HEj. Within mean 
field theory, it is possible to consider the effect of a non- 
zero temperature. For J±/J = 10 the temperature effect 
is negligible as long as T < 0.5J_l. For larger T, there is 
a decrease of s that causes a decrease of the bandwidth of 
triplet excitations (see Fig. ||). The effect of C(9i,9s) is 
to give an explicit dependence on the polarization of in- 
coming and outcoming radiation to the Raman intensity. 
However, it does not give rise to a shift of the peaks when 
the polarization is changed. In the case 7j_/7[| = J±/J\\, 
the Raman intensity is proportional to cos 2 (9] + 9s) as 
pointed out in Ref. |34|. 

Eq. (|24|) was derived neglecting triplet-triplet scatter- 
ing. Such scattering is r espon sible for the broadening of 
the peaks in experimentsGJii3't3ll3. A full calculation of 
the effect of triplet-triplet scattering is a very complex 
task. Physically, the most important effect of these in- 
teractions is that triplet bosons are not anymore exact 
eigenstates of the Hamiltonian and acquire a finite life- 
time r. An estimate of this lifetime is given by the Fermi 
Golden rule as: 



- = 2nJ2 \V{q)\ 2 n B {k'){\ + n B {k> + q)) 



k' ,q 



x S(e{k + q) + e(k -q)- e{k) ~ e(k')). 



(26) 



The potential for triplet-triplet interactions being 
short ranged, V(q) ~ V(q = 0) = Vq. Since the quasi- 
particle dispersion is of the form: ui k = uj-k + at the 
bottom of the spectrum, — will be dominated by a factor 
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r/T 



The rest of the expression reduces to a phase 



space factor, independent of temperature. 

The existence of a finite lifetime for the quasiparticles 
leads to the replacement 



S(uj - 2uj k ) — > - 



1 



7r( w -2 Wfe )) 2 + r 2 ' 



(27) 



in Eq. gjjj. At temperatures T -C c^o, T is small and 
such replacement does not affect the intensity very much 
except at the edges where it cuts the square root singu- 
larities of the density of states. Thus, we can restrict to 
the consideration of Raman intensity close to the edges 
of the spectrum. 

Let us first consider the threshold The Raman 

intensity is approximately: 



interesting to compare the dependence predicted by Eq. 
( |32j ) with the two magnon intensity of Fig. 4 in [l5| or 
Fig. 7 of Ref. [r| At temperatures small in comparison 
to the spin gap, the damping effect is the dominant fea- 
ture near the edges. At low temperature there is a rapid 
variation of the damping rate r = 1/r ~ e~ As/T with 
temperature and thus a rapid variation of the intensity. 
It would be interesting to check whether such dependence 
of the damping rate is consistent with experiments. As a 
final note, we would like to point out that Equations ( |30| ) 
and ( |32| ) do not lead in general to a symmetric Raman 
spectra. In particular, in the case of CaA^Os, a different 
damping rate and different masses of excitations could 
explain the. difference of behavior of the intensity at the 
two edges£3. 



I(u>) = coth 



coth 



\ 


( Ar 


k B T) 




LJir \ 


fK 


k B Tj 





dk r„ 



where 



/??„ 



TV J J 27T 7T (u) - 2u}ir 

2VT — 5 



k 2 \2 



. 4 



The calculation leads to: 



Six 2 



I(uj) = i Iid{u) = 



^x 2 + WVx 2 + l-x 



I (id) = coth 



where: 



mo ( 2cjo — u 

-IlD 



\2k B Tj \lj q J V 8^ 2 r LU \ r 



•n 



(28) 



III. COUPLED LADDERS 



A. Mean field equations 



(29) 



y/(LJ - 2A) 2 + Tl y/y/(u> - 2A) 2 + VI - cj + 2A 

(30) 

The resulting plot of I(u)) versus u> is shown on figure |^ 
and is qualitatively similar to the behavior of the Raman 
intensity close to 2A S in experimental systems (see Fig. 
3 in Ref. || or Fig. 5 in Ref. ||). Eq. © can be cast 
in the scaling form: 



> (31) 



where x = "p 2A . The corresponding plot is Fig. ||. It 
would be interesting to determine whether the experi- 
mental data satisfy such scaling form. A similar calcula- 
tion can be done for u> ~ ujq, with result: 



In this section we start to consider the effects on the 
Raman intensity due to interladder interaction. Our final 
aim is to derive an expression of the Raman intensity in 
the case of the array of ladders whose Hamiltonian was 
given m Eq. (|). Such ladder .-array model is of intec- 
est for systems such as KCuChjcS or possibly CuHpClcJ 
according to some neutron scattering experiments. As a 
simpler problem, we will also discuss a double rung ladder 
in which two spin- 1/2 ladders are connected by interlad- 
der exchange of strength X'J±, that is weaker than the 
intraladder coupling. In the BOT the Hamiltonian given 
by Eq. [li] where the index j can run now from 1 to N. 
First we apply the Green's function method to determine 
the spectrum, whose complete calculation is reported in 
the Appendix . It will be given by 

cu 2 (k x ,ky) = (-^- - \i - J\\s 2 cos(k x ) + - J ^ S cos(fcy) 



J||S 2 cos(k x ) — 



■ cos(ky) 



(34) 



The self-consistent equations in the saddle-point ap- 
proximation that permit us to determine the parameters 
(s 2 ,/u) at T=0 are given by 



(s 2 - -) 



71/(1 dk r r /a dk 



rr/a ( 27r ) J -ix /a (2tt) 

( 32 l3 ^ A r /a ^ r /a ^ 

[ 2 J± J-*/ a (2*r) J^ /a (2tt) 



1 



f(k x , ky) 



y/1 - df(k x ,k y ) 

f(k x . k y ) 
y/l-df{kx,ky) 



too 



2VT + d 



where 



(33) 



d = 



Although in experiments on Sri4Cu2404i the condi- 
tion J± ^> J is certainly not satisfied, it is nevertheless 



,2X 



f = (^ COSfc;,; - COSky) 



(37) 
(38) 
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The equations ( J35| ) and ( Pq ) can be combined to obtain 
the following single mean-field equation that is given by 



A. 



J± 
4 



l-d 1 



2 A 



1/2 



d 2 fc 



1 



(2tt) 2 77f • 



Once we have determined d, s 2 and the chemical po- 
tential /i, we are ready to perform the analysis of Raman 
susceptibility. 

In the case of the double ladder, the energy dispersion 
is determined by restricting k y to take the discrete values 
or 7T. This leads to two branches in the spectrum due to 
the fact that as soon as A' is different from zero the triplet 
excitations not only delocalize along the single ladders 
but also across the ladders. The two branches are usually 
called the bonding (the one below) and antibonding (the 
one above) band. 

The resulting excitation spectrum is given by 



<(/.•)= (^-M+ J|| S 2 cos(fc)±^ 



2N 2 



- ( J||S 2 cos(fc) =F 



J ± s 2 



(39) 



The chemical potential /i and the parameter s must 
obtained as above by solving the saddle-point equations, 
with the integral over k y changed to a discrete summa- 
tion. They read: 



(s 2 % + 



,3 „ u 
- + 2 — 

V 2 Ji 



where 



w ' a dk 



.*/„ (2tt) 
T/a _dk_ 

■tt /a ( 2?r ) 



As 2 



/+ 



f- 



= o 



It was shown by Gopalan et all 2 ! that as A' was increased, 
the spin gap was reducing. 

We can now turn to the calculation of the Raman sus- 
ceptibility. 

B. Derivation of the Raman Susceptibility 

For the calculation of the Raman susceptibility for cou- 
pled ladders we make use of the expression ( |l6|) for the 
Raman operator. To start with, we find that: 

(T t H r (t)H r (0)) 
= 3^ |A 2 (k) [G^G^-r) + F k (r)(f t) k (- T )] 

k 

+2A(k)B(k)cos(k x a) [G k (r)(Ft) k (-r) + G k (-r)(^) k (r) 

+ G k (-r) J F k (r) + G k (r)F k (-r)] 
+B 2 (k) [e 2 ^ Q (Ft) k (~r)(Ft) k (-r) +G k (r)G k (r) 

+ e - 2ife ^ k (T)F k (r) + G k (-T)G k (-T)] } . (45) 

Where G, F represent the normal and anomalous Green's 
Functions of the triplet bosons and are defined in the Ap- 
pendix. The term in A 2 (k) has already been calculated 
in the single ladder problem, while the remaining terms 
give new contributions to the Raman intensity. Going to 
Fourier space and making use of the Matsubara summa- 
tions reported in the Appendix, we find that the Raman 

response function is: 
(40) V 



o * R 



(iw n ) =EjE 2 s J2 M ^ coth 



1 



1 



2u k 



(*-£) 

/ ± = (cOS(fc)±™) 



(41) 
where: 

(42) M k = 4 J B^cos 2 (fc x a) 



(46) 



Ak 

w k 



2^1 k -B k cos(k x a) 



A k Ai 



(43) 



-2A 2 I — 



A k 
w k 



(47) 



Eqs. @ and @ can be combinedB to obtain the fol- 
lowing single mean-field equation 



In particular, if J\\/J\\ = Y_i_/ 'J± = l'x/J'x> °- ue ^° re l a ~ 
tion the expression for A/ k simplifies: 



d_ 2 2 1 

A 7T ^1 + ^(1 + a)' 



2d 



1 



K 



7T ^1 + d(l - a) \l + d(l-a) 



l + d(l + a) 
2d 



M k 



42 



A k Ai 




(48) 



(44) 



1 V 

2 A 



where a 

At this point it is worth analyzing the spin-triplet gap 
A s . It is given by 



where we have used the following notation: 



H) - I Jys COS^) 



cos(k y 



(2A k ) = J||S 2 cos(fc x ) 



A' J 1 



■ cos(fc !/ ). 



(49) 
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Note that the two coupled ladders case is included in the 
formula above since we can restrict the sum over k y to 
and 7r. The expression ( fl6| ) is very general and contains 
an explicit dependence on the polarization and the elec- 
tric field intensity. We emphasize that the peak position 
predicted by Eq. (^) only depends on the spectrum of 
the ladder system. In particular, it is completely inde- 
pendent of the polarization of ingoing as well as outgoing 
radiation. 

Performing an analytic continuation, the intensity of 
the Raman spectrum is given by 



d 2 k 

—-M k {5(LJ - 2cj k ) - 5{uj + 2w fcl )), 
(2tt) 2 



(50) 



where Mk is given by (j48|). Since the 7's appearing in the 
expression for Mk can be treated as phenomenological 
parameters we will perform the explicit calculation for 



iw/J\\ = i±/J± = iJA- 



1. The Raman intensity of the two-rung ladder 

Using the previous formulas, the Raman intensity for 
the two-rung ladder will be obtained from formula ( p30| ) 
reducing the integral on k y to a sum for k y = and 



Al(A k A k -Al) 



5(uj + 2u k , a ))- 



(5(w - 2oj k , a ) 
(51) 



Here we have summed on the two triplet bands. Per- 
forming the integral in (|5l|), we obtain as a final result 



(2(71/4-/.)) 



2 



2(7x74-/1) 



a=± 



4ui\ 



(jffe) (2(7774^1) a ° 



where An 



i = (J±/4-^) 3 (j) cos 2 (e I + e s ) ,a± = (l± 

As shown in Fig.^, the Raman spectrum shows four 
peaks in correspondence of the bottom and the top the 
bonding and antibonding bands. From simple density of 
states argument, no signal should be seen for lu < 2A S , 
where A s is the singlet-triplet gap. 



2. The array of ladders 

In this case we start from equation (|50|), by introducing 
the two-dimensional density of states p{oj) and we can 
write 



ImxflMoc ~p(oj/2)[( - — j -1], 



(53) 



/>M = -[77 



(^-m) J \[ 



2tt 



\ 



(54) 




2A COS 



The last integral is performed for Max(— 



2A 
dX' 



2A 



1 




< cosk v < MinCl, 



2A 
' A' 



The results show that the in- 

tensity is non-zero only between the maximum and the 
minimum of w(k) and a discontinuity appears at the band 
edge. The behavior of the 2D Raman intensity is sketched 
on Fig. ||. 

In two dimensions, we can expect two types of sin- 
gularities in p(i-lA, namely discontinuities and logarith- 
mic singularities^! Close to these singularities, damping 
effects due to triplet-triplet scattering will play an im- 
portant role, in analogy to the single ladder problem. 
In the present case the two-magnon continuum starts 

at energy: = 2A S = 2(J x /4 - n)yjl - d(l + Jo, 
the other edge of the Raman spectrum being at uj + = 

2(Jj_/4 — + d(l + 2^7). No magnetic Raman scat- 

tering is observed for oj ^ [u>-,ui+]. At these points, dis- 
continuities in the Raman intensity would appear in the 
absence of damping. Taking damping into account close 
to the edges of the spectrum, we find that the intensity 
is of the form: 



xodgc 
'2D 



, (52) 



kdk 



2ir (u - 2A - k 2 ) 2 + T 2 



arctan ■ 



2A 



(55) 



The steps at the edges of the spectrum are thus smoothed 
out by damping of excitations. The corresponding plot 
is Fig. I 

At a logarithmic singularity u> ~ 2(Jj_/4 — 
jj)Jl ± d(l - 2^7), we find that: 



xpeak 
L 2D 



(2tt) 2 K{uj-2A-k 2 x /m + k 2 y /m) 2 +T 2 ' 

(56) 



and thus: 



I{uj) ~ In 



A 2 



r 2 + (^-2A) 2 ' 



(57) 
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where A is an energy cutoff. Again, the height of the 
peaks is limited by the damping of excitations. Note 
that in contrast to the one dimensional case, the peaks 
are not obtained at the edges of the two magnon Raman 
spectrum but are superposed onto a positive two-magnon 
background. Only in the one dimensional limit do the 
peaks merge with the steps at the edge of the spectrum. 
Also, in the two dimensional case, the peaks are symmet- 
ric (see Fig. |l^) in contrast to the ID case. Thus, Raman 
scattering can prove useful to discriminate quasi-one di- 
mensional systeins-Jjrom quasi two-dimensional systems 
with a spin gap23'E2l and provide bounds on interladder 
coupling. 



IV. CONCLUSION 

We have presented an analysis of the Raman spec- 
tra for various coupled spin ladder systems based on 
the Loudon-Fleury photon-induced super-exchange the- 
ory, representing the spins using the Bond Operator For- 
malism. For a single ladder, we have shown that peaks 
should appear at the edges of the -Raman spectrum as has 
been observed in experimentstJ'Ej. The shape of these 
peaks is determined by the damping of triplet excitations 
due to interactions. A scaling plot of the Raman intensity 
versus frequency should permit to determine the lifetime 
of excitations in the ladder. Simple arguments lead us 
to expect that the lifetime of excitations must vary as 
e A / T where A is the spin gap. Another consequence of a 
positive temperature on the spin ladder is the narrowing 
of the spectrum as temperature is increased. It would 
be interesting to compare our predictions to the avail- 
able experimental data on spin ladder systems. In the 
case of the coupled array of ladders, we have shown that 
the Raman scattering intensity measures the density of 
states of triplet excitations. Thus, the singularities of the 
Raman intensity should indicate very clearly whether the 
system is one or two dimensional. In both case, the peak 
position was independent from the polarization of the in- 
coming and scattered radiation. It would be interesting 
to generalize our analysis to other two dimensional spin, 
gap systems such as the models with orthogonal dimersE£l 
that can be used to model SrCu2 (603)2- Another direc- 
tion in which-our work could be extended is to conducting 
spin laddersE3 in which case charge degrees of freedom 
could also come into play. 
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APPENDIX: GREEN'S FUNCTIONS METHOD 
FOR THE ARRAY OF LADDERS 

We introduce the following Green's functions: 

G Q (k,r)=-(T T t k , a (r)4 )a (0)) 
G Q (k,r) =-(T T ti a (r)t Ka (0)) 

F a (k,T) = -<T T tk,a(T)t_k,a(0)) 

(F^UKr) = -(T T tl a (r)tljO)) (Al) 

The equations of motion permit us to obtain the fol- 
lowing coupled equations: 

- iu! n G a (k,uj n ) = -1 - f^- - fi \ G a (k,uj n 

-s 2 cos(fc x ) - cos(k y )^ \G a (k, uj n ) + F^ (k, ui n )] 

-iuj n Fl(k,uj n ) = - [mJ i^(k,w„ 

+s 2 ^J|| cos(k x ) - ^y^cos(fc a )^ [G a {k,uj n ) + F+(k,o; n )] 



and: 

- iuj n G a (k, tj n ) = 1 + - f?j G a (k,cu ri 
+s 2 ^J|| cos(fc^) - cos(k y )j G a (k,cu n ) + F a (k,uj n ) 
-iuj n F a (k,uj n ) = - - f?j Fl(k, uj Ti 
-s 2 cos{k x ) - cos(fc,)) G a (k, u> n ) + F a (k,u n ) 



It is convenient to introduce: 

J J' s 2 

A(k) = — /i + J||.s 2 cos(k x ) ^— cos(k y ) (A4) 



and: 



J' s 2 

2A(k) = J|| S 2 cos(k x ) - -±— cos(fc y ) (A5) 



so that: 



G(k,. n) = iGM = -^^ 
F(k,. n ) = FHk,. n ) = -^ w 



(A6) 
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Where we have: 



2 (k)=A 2 (k)-(2A(k)) 2 



(A7) 



The dispersion for the isolated ladder is obtained by tak- 
ing J' ± = in the preceding formulas. 

A convenient decomposition of the Green's functions 

is: 



where C(Qi,6s ) is g iven by Eq. fl25|). 

Using Eqs. ( A12 ) one obtains the expression Eq. §24 
for the Raman intensity. 



1 



F\k,LU n ) = 



A(k) 

W (k) 



14 

1 - 
1 



A(k)\ 



1 



A(k)\ 



iuj n + w(k) 
1 



(k) ) iuj n - w(k) 
1 



■ W (k) 



(A8) 



Using this decomposition, the Matsubara sums are re- 
duced to: 



1 



X! ,7 



1 



1 



(3 ^ i(0J n + ^n) ± W(k) MJ n ± (j(k) 

, ! coth(^) 
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FIG. 1. The two dimensional array of ladders. 




FIG. 2. The double ladder 



X = 0.1 



s 2 (T) 




FIG. 3. The singlet order parameter s 2 (T) and d 2 (T) as a 
function of the temperature (in units of J±) for A = 0.1. 
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FIG. 4. The Raman intensity for fixed J/J± =0.1 and 
increasing T/J±. For T/Jj_ < 0.1 the effect of temperature 
on the Raman intensity is not visible. For T/J± > 0.5, the 
peaks are moving toward each other as a consequence of the 
diminution of s and A 
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FIG. 5. A plot of I(uj) versus uj/A for different values of 
r/A. Increasing temperature amounts to increasing F and 
results in a reduction of peak intensity and the apparition of 
a small intensity below the gap 
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FIG. 6. The scaling form of the one dimensional intensity 
Eq. ph 
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FIG. 7. Raman intensity for the two coupled ladders, for 
J/Jx = 0.3, J'/J± =0.1. 




FIG. 8. The behavior of the density of states in a two di- 
mensional system, showing a discontinuity at threshold and a 
logarithmic peak inside the spectrum. 
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FIG. 9. The behavior of Raman intensity at the thresh- 
old of the two-magnon spectrum in a two dimensional ladder 
array taking into account the damping of triplet excitations. 
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FIG. 10. The behavior of peaks in the Raman intensity for 
a two dimensional array with damping of triplet excitation 
taken into account. Note the symmetry of the spectrum, to 
be contrasted with Fig. |^. 
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